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Abstract 

We  study  the  complexity  of  approximating  the  smallest  eigenvalue 
of  a  univariate  Sturm-Liouville  problem  on  a  quantum  computer.  This 
general  problem  includes  the  special  case  of  solving  a  one-dimensional 
Schrodinger  equation  with  a  given  potential  for  the  ground  state  en¬ 
ergy. 

The  Sturm-Liouville  problem  depends  on  a  function  q.  which,  in 
the  case  of  the  Schrodinger  equation,  can  be  identified  with  the  po¬ 
tential  function  V .  Recently  Papageorgiou  and  Wozniakowski  proved 
that  quantum  computers  achieve  an  exponential  reduction  in  the  num¬ 
ber  of  queries  over  the  number  needed  in  the  classical  worst-case  and 
randomized  settings  for  smooth  functions  q.  Their  method  uses  the 
(discretized)  unitary  propagator  and  arbitrary  powers  of  it  as  a  query 
(“power  queries”).  They  showed  that  the  Sturm-Liouville  equation  can 
be  solved  with  0( log(l/e))  power  queries,  while  the  number  of  queries 
in  the  worst-case  and  randomized  settings  on  a  classical  computer  is 
polynomial  in  1/e.  This  proves  that  a  quantum  computer  with  power 
queries  achieves  an  exponential  reduction  in  the  number  of  queries 
compared  to  a  classical  computer. 

In  this  paper  we  show  that  the  number  of  queries  in  Papageorgiou’s 
and  Wozniakowski’s  algorithm  is  asymptotically  optimal.  In  particu¬ 
lar  we  prove  a  matching  lower  bound  of  n(log(l/e))  power  queries, 
therefore  showing  that  @(log(l/e:))  power  queries  are  sufficient  and 
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necessary.  Our  proof  is  based  on  a  frequency  analysis  technique,  which 
examines  the  probability  distribution  of  the  final  state  of  a  quantum 
algorithm  and  the  dependence  of  its  Fourier  transform  on  the  input. 


1  Introduction 

This  paper  deals  with  the  solution  of  the  Sturm- Liouville  problem  on  a  quan¬ 
tum  computer.  Quantum  computers  have  shown  great  promise  in  solving 
problems  as  diverse  as  the  discrete  problems  of  searching  and  factoring  mm 
and  the  continuous  problems  including  integration,  path  integration,  and  ap¬ 
proximation  lEnanEiEiizi.  The  main  motivation  for  quantum  computing  is 
its  potential  to  solve  these  important  problems  efficiently.  Shor’s  algorithm 
achieves  an  exponential  speedup  over  any  known  classical  algorithm  for  fac¬ 
toring,  but  until  the  classical  complexity  of  factoring  is  proven,  the  exponen¬ 
tial  speedup  remains  a  conjecture.  The  quantum  algorithms  for  integration 
provide  provable  exponential  speedups  over  classical  worst- case  algorithms, 
but  only  polynomial  speedups  over  classical  randomized  algorithms. 

Recently  Papageorgiou  and  Wozniakowski  introduced  a  quantum  algo¬ 
rithm  for  the  Sturm-Liouville  problem  m  which  uses  the  quantum  phase 
estimation  algorithm.  They  showed  that  quantum  algorithms  with  power 
queries^  achieve  a  provable  exponential  reduction  in  the  number  of  power 
queries  over  the  number  of  queries  needed  in  the  classical  worst-case  or  ran¬ 
domized  setting.  Naturally  query  complexity  results  neglect  the  cost  of  ac¬ 
tually  implementing  the  queries.  At  the  end  of  this  paper  we  will  discuss 
this  problem  for  power  queries,  but  it  is  currently  not  clear  under  which 
conditions  power  queries  are  sufficiently  inexpensive  to  implement  for  the 
Sturm-Liouville  problem. 

In  this  paper  we  will  prove  lower  bounds  on  the  number  of  power  queries 
for  quantum  algorithms  that  solve  the  Sturm-Liouville  problem.  This  can 
be  used  to  show  the  optimality  of  the  algorithm  proposed  in  hi.  To  prove 
lower  bounds  for  algorithms  with  power  queries  the  previously  known  quan¬ 
tum  lower  bound  techniques,  such  as  the  “polynomial  method”  of  Beals  et.  al 
mini  do  not  suffice.  Our  lower  bound  method  builds  on  the  “trigonometric 
polynomial  method”  [2],  which  is  an  extension  of  the  above-mentioned  poly¬ 
nomial  method  and  was  modified  to  be  used  with  power  queries  in  [3]  to  prove 

1We  will  define  power  queries  rigorously  in  Definition  [l]  Informally  they  are  just  an 
arbitrary  (integer)  power  of  a  specific  unitary  matrix. 
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lower  bounds  for  the  phase  estimation  algorithm.  Our  method  uses  frequency 
analysis  instead  of  a  maximum  degree  argument,  which  is  not  applicable  in 
the  case  of  arbitrary  powers. 

2  The  Sturm-Liouville  eigenvalue  problem 

Papageorgiou  and  Wozniakowski  study  in  HI  a  simplified  version  of  the 
univariate  Sturm-Liouville  problem.  Consider  the  eigenvalue  problem  for  the 
differential  equation 

—u"{x)  +  q(x)u{x)  =  A  u{x) 

i/(0)  =u(  1)  =  0  ^  ) 

for  a  given  nonnegative  function  q  belonging  to  the  class  Q  defined  as 
Q  =  {q  ■  [0, 1]  -*•  [0, 1]  :  qe  C2([0, 1])  and  max  max  |gW(x)|  <  l}.  (2) 

We  are  looking  for  the  smallest  eigenvalue  A  such  that  there  exists  a  non¬ 
zero  function  u\  that  satisfies  ([Tj) .  What  is  the  minimal  number  of  queries  of 
q  that  permits  the  determination  of  the  smallest  eigenvalue  A  in  this  equation 
with  error  e  and  probability  3/4  on  a  classical  or  quantum  computer? 

The  one-dimensional  time-independent  Schrodinger  equation 

h2  d2 

- —  Tbr)  +  Wz)T(x)  =  EV(x)  (3) 

2m  ax2 

of  a  particle  in  a  box,  see  Pi.  is  an  instance  of  ([Tj) .  We  are  given  a  potential 
V  and  are  looking  for  the  eigenfunctions  T  of  this  equation  and  their  corre¬ 
sponding  energies  E.  In  particular,  we  are  interested  in  the  ground-state  and 
its  energy,  i.e.,  for  a  given  potential  V,  we  want  to  determine  the  eigenfunc¬ 
tion  To  and  its  energy  Eq,  such  that  all  other  eigenfunctions  Tn  have  higher 
energies  En  >  E0.  Since  quantum  systems  obey  equation  (|3]),  it  seems  plau¬ 
sible  that  quantum  computers  could  potentially  solve  the  eigenvalue  problem 
faster  than  a  classical  computer. 

In  the  next  section  we  define  a  quantum  algorithm  with  power  queries. 
We  especially  have  to  tackle  the  question  concerning  the  form  of  the  input 
(i.e.,  the  function  q  in  the  Sturm-Liouville  problem)  enters  the  quantum 
algorithm. 
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3  Quantum  algorithms  for  the  Sturm-Liouville 
problem 


Let  us  denote  the  differential  operator  associated  with  the  Sturm-Liouville 
problem  for  a  certain  q  G  Q  as  hq  :  C2([0,  1])  — ■>  C°([0,  1]),  defined  by 

d2 

L  „u(x)  =  — — 7u(x)  +  q(x)u(x). 
axz 

We  discretize  L q  by  approximating  the  second  derivative  at  the  points 
— tt ,  . . -f-r  and  obtain  an  n  x  n  matrix  Mn: 

n+1  ’  ’  n+1  h 


'2-1 

-1  2  -1 

9(^l) 

Mq  —  (n  +  1)2 

-1  2  -1 

+ 

—  !  2  _ 

L  ?(^t)J 

The  eigenvalues  of  L9  and  Mq  are  closely  related.  Let  us  denote  the  smallest 
eigenvalue  of  L„  by  X(q)  and  let  us  write  Ai  (Ma)  for  the  smallest  eigenvalue 
of  Mq.  Then  (see  e.g.  0) 

\(q)  -  X^Mq)  =  0(n~2).  (5) 


The  input  q  G  Q  enters  the  quantum  computer  in  the  form  of  a  unitary 
black-box  transformation  called  a  quantum  query.  For  the  Sturm-Liouville 
problem  we  define  this  query  to  be  the  unitary  operator  exp (|Mg).  One 
can  show  that  the  smallest  eigenvalue  A (q)  of  the  Sturm-Liouville  equation 
satisfies  7r2  <  A (q)  <  7r2  +  1.  To  avoid  ambiguity  we  use  proper  scaling,  i.e. , 
instead  of  exp (iMq)  we  use  exp(|M9),  which  dehnes  a  unique  phase  ip  G  [0,  1) 
by  2ni<f  =  |A (q). 

We  now  dehne  an  associated  quantum  power  query  for  exp (|Mg). 

Definition  1.  Let  hq  be  the  differential  operator  for  a  Sturm-Liouville  prob¬ 
lem  and  Mq  its  discretization  at  n  points  as  in  |^)  for  q  G  Q.  We  define  the 
power  query  Wf  (exp(^ Mq))  as 


Wf(exp(|Mq))  |xi) . . .  \xc)  \  fj) 


|aq)  . . .  \xc)  exp {\pMq)  \fj)  for  xt  =  1 
|xi) . . .  \xc)  | if)  otherwise 


for  all  xi,...,xc  G  {0,1}  and  arbitrary  normalized  vectors  \if)  G  Cn  and 
extend  this  definition  to  all  quantum  states  by  linearity. 
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Suppose  that  the  \ifs),  s  =  1, . . .  ,n,  are  the  eigenvectors  of  Mq  and  that 
Mq  | ips)  =  \s  \ips).  Then  for  \ip)  =  Y^=ias  \i>s)  and  |x)  =  |aq) . . .  \xc)  with 
xi  =  1 

n 

Vhf(exp(|Mg))  |x)  \if)  =  |x)exp {\pMq)  |V>)  =  |x)  e2pXs  \ips) . 

s= 1 

Quantum  algorithms  are  products  of  unitary  transformations.  Every 
quantum  algorithm  that  approximates  A  (q)  can  be  divided  into  stages  that 
use  powers  of  exp(|Mg)  and  therefore  depend  on  q,  and  stages  that  are  in¬ 
dependent  of  q.  Let  us  define  a  quantum  algorithm  with  power  queries. 


Definition  2.  For  a  Sturm- Liouville  problem  given  by  the  input  q  G  Q  with 
the  solution  A  (q),  we  define  a  quantum  algorithm 

A  =  Uq,  ■  ■  ■ ,  Ut ;  h,Pi,  ■  ■  ■ ,  It,Pt,  A) 


with  T  power  queries  that  solves  this  problem  as  follows.  Let  U0,  U\,  . . Ut 
be  arbitrary  but  fixed  unitary  transformations  and  a  fixed  initial  state. 


Let  Wf3  (exp(|Mq))  be  a  power  query  as  in  Definition 
the  state 


A  measurement  of 


p(T>(exp(|M,)))  =  UTWp(exV(iMq)) . . .  lAWf? (exp(iM,))C/„|V>(<,)> 

in  the  standard  basis  yields  a  state  |  k)  with  probability  Pk(q)-  For  each  k 
compute  an  approximation  A (k)  G  M  to  the  eigenvalue  of  interest  A (q)  on  a 
classical  computer.  For  every  q  G  Q  the  probability  that  an  e- approximation 
A (k)  of  X(q)  is  computed  is  given  by 

J2  Pk (?)■  (6) 

fc:|A(<j)  —  A(fc)|<£ 

For  any  algorithm  A  with  T  power  queries  we  define 

e(A,T)  =  inf  {e  :  e  chosen  such  that  (jb|j  is  larger  than  |  for  all  q  G  Q} 

as  the  worst-case  quantum  error  of  A. 

We  measure  in  the  standard  basis  for  convenience  only;  a  measurement  in 
any  other  basis  is  easily  achieved  by  modifying  the  operator  Ut  accordingly. 

A  model  like  this  was  introduced  in  [lj  for  discrete  inputs  q.  It  was  ex¬ 
tended  to  continuous  functions  by  Heinrich  in  [5].  Our  model  is  an  extension 
of  this  model  to  incorporate  power  queries. 
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4  Upper  bounds 


To  estimate  A (q)  on  a  quantum  computer  with  power  queries  Papageorgiou 
and  Wozniakowski  used  the  quantum  phase  estimation  algorithm,  see  e.g. 
[12].  This  algorithm  takes  a  unitary  transformation  Q  with  an  eigenvector 
|£)  as  input,  i.e. ,  Q  |£)  =  e2mtp  |£).  Here  p  G  [0,  1)  is  called  the  “phase”  of  the 
eigenvalue  corresponding  to  |£),  and  the  phase  estimation  algorithm  gives  us 
an  approximation  ip  to  p.  This  algorithm  has  the  final  state 


|V>(T)«))  =  (Tp  ®  i)wr\Q)wr\Q)  ■  ■  ■  w?(Q)(h®t  0  /)  |o>  |0 , 

and  is  depicted  in  Figure  [l]  Suppose  Q  is  a  r  qubit  transformation.  A 


|0> 

|0) 

|0) 

|0) 

10 


10) 

10) 

|^T-l) 

\kT) 

10 


Figure  1:  The  quantum  phase  estimation  algorithm.  T.2t  is  the  inverse 
quantum  Fourier  transform  on  T  qubits. 


measurement  of  \^T\Q))  returns  a  state 

I  k)  =  \h) . . .  I  hr)  \kT+i)  ■  ■  ■  |  kT+r)  ■ 

The  algorithm  then  uses  k  to  compute  an  approximation  p{k)  =  k\2~x  + 
k22~2  +  . . .  +  kx 2~T  to  p  classically. 

One  can  show,  see  e.g.  [12|,  that  with  probability  greater  than  |  the  al¬ 
gorithm  approximates  p  up  to  precision  e  with  C,(log((l/e)))  power  queries. 
Papageorgiou  and  Wozniakowski  use  this  algorithm  to  approximate  the  small¬ 
est  eigenvalue  A (q)  of  the  Sturm-Liouville  operator  hq  and  use  the  operator 
Q  =  exp(^Mg)  as  a  query.  Since  the  phases  of  exp (|M?)  and  exp(|L?)  are 
related  through  equation  ([H]) ,  we  have  to  discretize  at  n  =  (P(£“1,/2)  points. 
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The  quantum  phase  estimation  algorithm  requires  the  knowledge  of  the 
eigenvector  for  which  the  phase  is  estimated.  For  the  Sturm-Liouville  prob¬ 
lem  we  need  the  eigenvector  \z\(Mq))  of  Mq  corresponding  to  the  small¬ 
est  eigenvalue  Xi(Mq).  We  can  compute  \z\{Mq ))  through  the  method  of 
Jaksch  and  Papageorgiou  [8],  which  computes  a  superposition  of  eigenvec¬ 
tors  | Zj(Mq))  of  Mq ,  with  a  sufficiently  large  \zi(Mq))  component.  For  details 
see  (HI  LB]  • 

5  Lower  Bounds 

Our  goal  is  to  prove  that  the  algorithm  described  in  the  previous  section 
is  optimal  with  respect  to  the  number  of  power  queries.  We  have  to  prove 
that  every  quantum  algorithm  A  with  T  power  queries  that  returns  a  correct 
answer  with  precision  e(A,  T)  <  £  has  to  use  T  =  fi(log(l/e))  power  queries. 

We  will  show  that  even  for  a  much  simplified  version  of  the  problem  this 
lower  bound  still  holds.  Consider  as  input  only  constant  functions  q{x)  =  q  G 
[0,  1].  Obviously  q  G  Q.  It  is  easy  to  see  that  in  this  case  the  eigenfunctions 
which  fulfill  the  boundary  condition  in  ([Tj)  are 

us(x)  =  sin(s7nr)  (7) 

for  s  G  N  and  that  they  have  eigenvalues  As  =  s2tt2  +  q,  which  means  that 
the  smallest  eigenvalue  A (q)  is  A (q)  =  tt2  +  q. 

Similarly  for  the  discretization  Mq  of  hq  with  constant  q  G  [0,  1]  the 
eigenvectors  are 

n 

\u»)  =  \/4i5Zsin(iwr)  l*>  (8) 

x=l 

with  eigenvalues  4(n  +  l)2  sin2  (yyyy,)  +  <1- 

We  want  to  investigate  how  different  power  queries  lead  to  different  out¬ 
puts  and  turn  to  the  techniques  in  |3|.  For  the  following  theorem  we  abbre¬ 
viate  j  =  (j i, . . .  ,jn),  if  =  (cpi, . . . ,  (pn),  and  use  the  following  “dot  product” 
notation: 

J  ■  <P=jl<Pl  +  ■  ■  ■  +  jn<Pn- 

Theorem  3.  For  a  given  orthonormal  set  of  vectors  {| ifs)  :  s  —  1,  2, . . . ,  n} 
consider  all  unitary  matrices  Q  for  which  the  |^s)  are  eigenvectors.  Denote 
the  eigenvalues  corresponding  to  \ ipf)  by  e2niip3 ,  ip8  G  [0,  1),  s  —  1,2, ...  ,n. 
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Any  quantum  algorithm  with  power  queries  W^Q),  see  Definition  ^  that 
uses  c£N  control  qubits,  can  be  written  as 


n  2C-1 


|V>(T)(Q)>  =  UtW%(Q)Ut-,  . . .  u,w™  W)C/„|^<0)>  =  £  S<T>(cp)  I k) , 


k= 0 


where  are  trigonometric  polynomials  of  the  following  form: 


(9) 


sPw  =  v  y™*,  (io) 


jeJr 


with  Jt  defined  as  J0  —  {(0, . . . ,  0)}  and 

Jt+1  —  {(ill  ■  •  •  j  Jn)  j  (il  +  Pt+  1)  ■  •  •  5  in)  5  ■  •  •  J 

(il,  •  •  •  ,jn+  pT+ 1)  :  (il,  •  •  •  ,jn)  e  Jt },  (11) 


and  the  coefficients  a'k  1  E  C  do  not  depend  on  ip  and  are  normalized: 


EE 

fc  JJt 


a 


CO  12  _ 


=  1. 


(12) 


The  proof  of  Theorem  [3]  can  be  found  in  the  appendix.  It  is  very  similar  to 
a  proof  in  [3],  but  additionally  proves  equation  (12). 

In  the  Sturm-Liouville  problem  we  are  interested  in  the  smallest  eigen¬ 
value  A (q)  of  Lf/.  Suppose  we  choose  the  constant  functions  q{x)  =  q  E  [0,  1) 
as  input.  The  eigenfunctions  us  of  hq  and  the  eigenvectors  |ns)  of  Mq  are 
then  independent  of  q,  see  ([7]),  ([8]).  Therefore  we  can  apply  Theorem  [3]  and 
write  the  final  state  | fi>^T\q))  =  |^i^(exp (|M?)))  of  any  algorithm  with  T 


power  queries  as 


|^(T)(<?)}  =  £5f(fl  \k)  =  £  £ 

k  k  j£jT 


a(T}e2ffi-0 

k,j 


\k )  i 


see  equations 
are  4(n  +  l)2  sin 
are 


(|9j)  and  (floj), 
^  (^rr£+ 


with  <p  depending  on  q.  The  eigenvalues  of  Mq 
q.  This  implies  that  the  eigenvalues  of  exp  (|M?) 


2 mvs  =  ei(4(n+1)2sin2  (  2(7+1)  )+«). 
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Therefore  we  can  split  e 2mLv  as  follows: 

e2 nij-0  =elG'1  (4("+1)2sin2  (  2(n+l)  )  +9)  +",+J"n  (4(«+l)2  sin2  (  2(^+1)  )  +<?)  ) 

=e|(ii4(n+1)2Sm2  (  2(^1)  )+-+Jn4(n+l)2  sin2  (2(^1))  )  (h+-+in)<7 

_ A.  .  ,pf(jl+-+in)« 

which  greatly  simplifies  V^(9))  as 


|V-(t>(9)>  =E 


a 


(T) 


(u, . we!<*+-+««  |fc) 


0'l>— Jn)eJr 

-EE 

fc  itiGMt 


(13) 


with 

•M.T  —  {jl  +  •  •  •  +  in  :  (jl>  •  •  •  T.'in)  €  Jr}  (14) 

and  the  coefficients  77^  defined  by 


Vk,m.  ~  E  ak,(j1,...,jrl)C(jl,-Jn)- 

jl+...+jn=m 


These  coefficients 


are  bounded  in  the  following  way: 


E  E  I7©2 

k  tti^AAt 


< 


EE  E  l4T(k,..hn)C0l,..hn) 


k  rn£MT  (jl,-,jn)€jT 

jl+...+jn=m 


EE 

k  J€Jt 


a 


cn.2 


=  1, 


(15) 


where  we  used  \C(ji,...,jn) 


1  and  equation  (12). 


Lemma  4.  Let  A  be  a  T  power  query  quantum  algorithm  for  the  Sturm- 
Liouville  problem  with  powers  p±, . . .  ,pr  and  c  G  N  control  bits  as  defined  in 
Definition  [£|  Let  B  be  a  partition  of  the  set  of  all  basis  vectors,  i.e. 


U  B  =  {\k)  :  k  =  0,1,...,  n2c-  1}  and  B  D  C  =  0  for  B,C  G  B,  B  ±  C. 

BeB 
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If  the  input  q  E  Q  is  a  constant  function  q{x)  —  q  E  [0,  1)  ,  the  probability  of 
measuring  a  state  \k)  from  B  E  B  is  a  trigonometric  polynomial 

Ml)  =  £  /3gM\  (16) 

1&CT 

(T) 

with  coefficients  (yB  /  E  C  that  are  bounded  by 

EiA3i<i 

BeB 

for  all  possible  partitions  B,  and  the  set  Ct  is  given  by  £0  —  {0}  and 

Bt+i  =  [J  {fl  +  Pt+iJ  ~  Pt+i}  ■  (17) 

i&cT 

Lemma  |4]  is  proven  in  the  appendix.  Note  that  \Ct\  <  3T.  This  bound  is 
sharp,  since  for  the  choice  of  Pi  =  3*_1  we  have  £0  =  {0} ,  C\  =  {  —  1,0,1}, 
£2  =  { — 4,  —3,  —2,  ...,3,4}  and  in  general 

CT  =  {— 3T  -  S7”1  -  . . .  -  1, . . . ,  3r  +  S7”1  +  . . .  +  1}  . 

5.1  Fourier  Analysis  of  Power  Query  Algorithms 

With  Theorem  [3]  and  Lemma  [4]  we  have  the  tools  needed  to  provide  a  lower 
bound  for  the  Sturm-Liouville  problem.  We  are  now  able  to  apply  our  fre¬ 
quency  analysis  technique  to  this  problem. 

Theorem  5.  Any  quantum  algorithm  A  with  T  power  queries  which  esti¬ 
mates  the  smallest  eigenvalue  A  (q)  in  the  Sturm-Liouville  eigenvalue  problem 
for  all  inputs  q(x)  =  q  E  [0,  1)  with  precision  e(A,  T)  <  e  and  probability 
greater  than  3/4  has  to  use  T  =  f2(log(l/e))  power  queries. 

Notice  that  a  lower  bound  on  the  “easy”  subset  of  constant  functions 
q(x)  =  q  implies  that  the  same  lower  bound  holds  for  any  set  of  inputs  that 
includes  the  constant  functions,  hence  it  also  holds  for  the  class  Q. 

Proof.  After  T  power  queries  we  measure  the  final  state  and  receive  a  state 
| k)  with  probability  Pk(q)-  From  the  integer  k  we  classically  compute  a 
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solution  A (k).  A  successful  algorithm  has  to  return  an  ^-approximation  for 
every  q  G  [0,  1)  with  probability 

Pk (?)  A 

fc:  |  A(qr) — A(fc)|<e 


see  Definition  [2}  Define 

Aq,e  :=  {k  :  |A(g)  -  \{h)\  <  e} 


as  the  set  of  states  that  are  mapped  to  ^-correct  answers  for  input  q.  Choose 
iVeN  such  that  A  is  slightly  bigger  than  2e,  i.e.,  <  2e  <  A  and  define 

the  points  xr  :=  (r  +  1/2 )/N  for  r  =  0, 1, . . . ,  N  —  1.  For  the  inputs  q  =  xr 
we  can  visualize  the  quantum  algorithm  A  as  in  Figure  [2j  Notice  that  the 
sets  AXr:£  are  mutually  disjoint  for  r  =  0, . . . ,  N  —  1,  because  xr  and  xr+\  are 
chosen  such  that 


|A(xr)  -  A(xr+i)| 


16  sin2  (|)  +  AA-  -  16  sin2  (f) 


r±i±j 

N 


1 

N 


>  2e. 


Therefore  there  can  be  no  state  |  k)  that  is  mapped  to  an  output  A (k),  which 
is  an  ^approximation  to  A(ay)  and  A(av+i)  at  the  same  time.  Let 


PrM)  =  Pk (9) 

k(zAxr,E 


(18) 


be  the  probability  of  measuring  an  ^-approximation  to  X(xr).  Since  the  sets 
AXr:E  partition  the  set  of  all  outputs,  Lemma  [4]  allows  us  to  write 

PrM)  =  PAXrM)  =  €]Mlq- 

1&CT 

We  apply  the  W-point  inverse  discrete  Fourier  Transform  to  priS{q),  which 
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Input  q 


Figure  2:  A  quantum  algorithm  for  the  Sturm- Liouville  problem  with  inputs 
q  =  xr,  r  =  0, . . . ,  N  —  1,  will  result  in  a  probability  distribution  Pk(q)  on  the 
states  |  k)  that  are  measured.  Each  state  |  k)  is  mapped  to  an  answer  A (k).  We 
write  AXr>e  for  the  set  of  all  states  |  k)  that  are  mapped  to  ^-approximations 
of  A(xr). 


we  evaluate  at  the  points  xn,  and  get  the  following  value  at  k  —  0, 1, . . . ,  N—l: 


N- 1 

DFTN\prjE\(k)  =Y,PrAxn)e-^ikn/N 

n= 0 
N-l 

—  p(T)^e\l{n+l/2)/N e~2mkn/N 

fi= 0  /g£t 


N-l 


E 

e2m(^-k)n/N 

1&CT 

n= 0 

■ 

f  e27ri(3V_fe)-l  l 

E  €W,(2N> 

iecT 

<  e2«(^-fc)/JV_i  ’  4tt 

1  A"  — 

t  J  4vr 

^  k  (mod  N ) 
=  k  (mod  N) 


(19) 


where  =  k  (mod  N)  indicates  that  there  exists  an  integer  z  such  that 
=  k  +  zN.  For  every  l  define  1/a-k{n)  G  [0,  N)  as 


1/a-k{n)  :=  min{Z/(47r)  —  zN  :  z  —  0,1,2,...,  and  l/(4n)  —  zN  >  0}  . 


Then  exp (271 n^/N) 
equation  (19),  we  use 


=  exp(27r U/att{n)/N).  To  take  the  absolute  value  of 
I el2e  —  1 1  =  2  |sin(0)|  and  get 


\DFTN\pr,£](k)\  < 


ieCT 


\3{t) 

\Pr,e,l 


|sin(7r(i/4T(^)  — fc))|  . 

|sin(7rP/47r(jV)-fc)/V)|  ’ 


N 


,1 


/Air(N) 


=  k 


(20) 


We  can  bound  the  Fourier  transform  (19)  by  separating  the  correct  an¬ 
swers,  i.e.,  the  e-approximations  to  ay,  from  the  rest:  if  the  input  q  =  ay 
then  the  algorithm  has  to  return  an  answer  A  that  is  e-close  to  the  correct 
answer  A  (ay)  with  probability  greater  or  equal  3/4.  This  probability  is  given 
by  pr,e(q)i  be.,  we  demand  that  pr,£(ay)  >3/4.  Then: 


N- 1 

^  Pr,e(Xn)e 

n= 0 


N—l 


—2irikn/N 


> 


Pr,e(Xr )  ^  ^  Pr,e(^n) 


n=0 

n^p-r 


3 

>- 

~4 


N-l 

'FVt~ 


(21) 


Xn)i 


n= 0 
n^r 


Consider  the  second  term  in  (21),  ^2n=o  Pr,e(xn)-  Recall  that  pr,£(q )  is  the 

n^r 

probability  that  the  algorithm  measures  a  state  |  k)  that  is  mapped  to  an 
answer  A (k)  that  is  an  e-approximation  to  A(ay),  i.e.,  | k)  G  AXr>£,  see  (18). 
This  probability  pr,e(q)  depends  on  the  actual  input  q.  For  input  q  =  xn  ^  ay 
a  state  | k)  G  AXri£  will  not  yield  an  e-correct  answer:  we  chose  the  xn, 
n  —  0, . . . ,  N  —  1,  such  that  |A(ay)  —  A(ay)|  >  2e  for  n  ^  r,  and  thus  there 
cannot  be  an  e  close  answer  for  both  xr  and  xn.  The  sum  Y/n=(i  Pr,e(xn )  now 

n^r 

tells  us  how  often  the  algorithm  chooses  a  state  from  Are. 

If  we  knew  that  none  of  the  wrong  answers  is  preferred  by  our  algorithm, 
say  e.g.  J2n=o  Pr,e(xn)  <  §,  equation  (21)  would  read 

n^r 


N-l 

y  PrF 


Xr, 


n= 0 


3 

>  - 
-  4 


N-l 

yi  Pr,e  I 


Xr, 


n= 0 
n^r 


l 

>  4 


(22) 
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We  will  show  that  this  property  has  to  be  true  for  some  r  =  0, . . . ,  N  —  1, 
indexing  the  set  of  states  AXrj£  that  represents  numbers  e-close  to  xr.  Let 
R<  be  the  set  of  all  r  for  which  J2n=o  Pr,s(xn)  <  \  holds  and  R-  the  set  for 

riy^r 

which  it  does  not.  We  estimate  the  number  of  elements  of  i?<  by  splitting 

N- 1  N-1N-1  N-l  N-1N-1 

R  ~  ^  ^  1  A  y  }  Pr,e{xn)  =  ^  ]  Pr,e{xr)  +  ^  ^  ^  ^  Pr,e{xn) 

72—0  72—0  r=0  r=0  r=0  n=0 

n^r 

into  the  following  parts: 


N-l 


N-l 


N-l 


N  R  y  }  Pr,e(^r)  A  ^  ^  }  Pr,e(xn)  A  ^  ^  ^  Pr,e{p^n) 


r= 0 


reR<  n= 0 
n^r 


+  |i?<|  -o 


r£R>  n= 0 
n^r 


and  therefore  we  can  conclude  that  | R- 1  <  \N  and  thus  |i?<|  >  ^N.  Now 
\R<\  >  0  implies  that  we  can  actually  choose  an  element  r  G  R<.  Fix  such 


an  r  and  we  can  combine  equations  (20)  and  (22)  to 


V4  <  £  |#; 


(T) 

£,l 


|sin(7r(;/47r(Ar)-fc))| 

|sin(7r(Z/4w(jv)— fe)/AT)  |  ’  /AAN) 

A"  ,  l/4n(N)  k 


7^  k 


(23) 


We  will  now  fix  the  parameter  k  =  0, 1, . . . ,  N  —  1  in  inequality  (23)  in  such 


a  way  that  the  terms  in  the  sum  on  the  right- hand-side  of  the  inequality 
are  as  small  as  possible.  This  will  imply  that  the  sum  must  be  over  a  large 
number  of  elements,  i.e. ,  that  \Ct\  is  large.  Since  \Ct\  <  3T  this  will  help  us 
to  ultimately  prove  that  T  =  fi(logiV)  if  we  could  show  that  \Ct\  —  0(1V). 
More  specifically  we  will  show  that  \£t\2  >  ttW  which  proves  T  =  fl(logiV). 


We  prove  \Ct\  >  jqN  by  contradiction.  Assume  \Cq 


< 


A  TV 
10  • 


This 


assumption  allows  us  to  hnd  a  k  such  that  the  right-hand-side  of  inequal¬ 


ity  (23)  is  smaller  than  the  left- hand-side,  which  will  lead  to  our  desired 


contradiction. 

If  we  project  Ct  into  the  interval  [0,  N )  through  l  >  l /^(N)  we  will  get 
a  set  {//^(iv)  :  l  G  Ct}-  Order  this  set  as  0  <  tx  <  t2  <  . . .  <  t\CT\  <  N. 
This  defines  “gaps”  between  these  numbers,  i.e.,  intervals  G  =  (tj,  tj+ 1)  for 
j  =  1, . . . ,  \Ct\  if  we  dehne  t|£T|+i  =  t\  +  N  (we  “wrap  around”).  Dehne  the 
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width  w(G)  of  such  a  gap  G  as  the  distance  between  its  endpoints.  Thus 
tj+ 1))  =  tj+ 1  -  tj. 

Let  Gm  be  the  gap  with  the  maximal  width  w(Gm )  in  the  distribution. 
Its  width  must  be  w(Gm )  >  N/  \Ct\ ,  since 


IV  =  max  in  (G)  =  |£t|  max  in  (G). 

G  G 

Additionally  in(Gm)  >  10,  since  we  assumed  \jCt\2  <  and  therefore 
i^i  >  10  \Ct\  >  10.  Thus  there  are  at  least  ten  integers  k  G  {0, 1, . . . ,  N  —  1} 
that  fall  into  this  largest  gap  Gm,  i.e  k  G  Gm.  One  of  these  k  has  maximum 
distance  to  both  boundaries  tj  and  tj+ 1  of  Gm:  it  is  the  k  that  is  closest  to 
the  middle  m  =  tj+i+L  of  Gm  =  ( tj ,  tj+i).  This  integer  k  fulfills  | k  —  m\  <  \ 
and 


min  {k  —  tj,  tj+i  —  k}  =  min  {m  —  tj  +  k  —  m,m  —  k  +  fJ+1  —  m} 

>  w{Gm)  _  1  >  IV  _  1 

~  2  2  “  2|£t|  2' 

Fix  this  k  G  (b,-,  tj+i)-  Now  |sin(x)|  >  2/7T  |x|  for  —  7t/2  <  x  <  7t/2  and 
therefore 


min 

/e£T 


sin 


/&ir(N)  k'j 

N 


>  min  —  1 1 
iecT  N 


/Att{N) 


1  1 


k 


—  min  {k-tj,tj+1  - 


k} 


Then  we  can  use  this  to  estimate 


1/4  <  E 

IGlCt 


|sin(7r(//47r(jv)  -  fc))| 

|sin(7r(I/47r(jv)  -  k)/N)\ 


sEh 

i&cT 


(T) 

r,e,l 


1/  \Ct\  -  1/N 


We  sum  the  last  inequality  over  all  r  G  R<  for  which  it  is  valid,  and  get: 


v  1  < - 1 - 

ir< 4  -  v  \ct\  -  i/jv 


E  E  |C 

rS-R<  1&CT 


(24) 


Since  the  number  of  elements  in  R<  is  bounded  by  |A<|  >  |IV,  the  left-hand- 
side  of  (24)  is  bounded  by  |IV  <  |I?<|  The  right-hand-side  of  inequality 
(24)  can  be  bounded  through  Lemma  111 


—  \£'t\ 


E  E  K 

rGi?<  IGlCt 
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1 2  1 

If  we  put  both  sides  together  again  and  recall  that  we  assumed  \CT\  < 
we  get 


1 

8 


N  < 


|  3 

Ct\ 

| 

Ct\ 

2 

1/  \£i 

■1/N  1  - 

\Ct 

\/N 

10 


N 


1  - 


i 

io|£t| 


< 


To* 

!-hi 


1 

9 


N, 


which  is  a  contradiction. 

Therefore  |£T|2  >  jqN  must  hold.  This,  together  with  \CT\  <  3T ,  leads 
us  to  N  <  10  •  9t.  Take  the  logarithm  and  we  get  T  =  fl(loglV).  We  chose 
N  such  that  <  2e  <  A  which  finally  proves  that  the  number  of  power 
queries  T  for  any  algorithm  A  with  error  e(A,  T)  <  e  has  to  be  of  the  order 
T  =  fi(log(l/e)).  □ 


6  Discussion 

In  this  paper  we  have  proven  lower  bounds  for  the  number  of  quantum  power 
queries  for  the  Sturm-Liouville  problem  and  settled  an  open  problem  in  da. 

How  does  this  number  of  T  =  ©(log(l/e))  power  queries  relate  to  the 
cost  of  quantum  algorithms?  Here  we  understand  “cost”  as  an  abstraction 
on  the  number  of  elementary  quantum  gates  or  the  duration  for  which  a 
Hamiltonian  has  to  be  applied  to  a  quantum  system.  Suppose  the  function 
q  is  from  a  class  Q'CQ  where  each  power  query  H//p(exp(|Mg))  can  be 
implemented  with  cost {Wf  (exp (|Mg)))  =  cost(Q',p). 

If  we  implement  IHp(exp(|M,?))  naively  as 

^(exp(|M,))  =  (HY(exp(|M,)))p, 

then  cost(Q/,p)  =  p  ■  cost(Q')  and  the  cost  of  the  Sturm-Liouville  algorithm 
with  T  =  @(log(l/ s))  power  queries  grows  as 

T—l  T—l 

cost(Q/,  2J)  =  •  cost(Q/)  =  (2r  —  1)  •  cost(Q')  =  ©(1/e)  •  cost(Q/). 

i= o  3=0 

This  is  polynomial  in  1/e  just  like  the  Sturm-Liouville  algorithm  with  bit 
queries  discussed  in  p3j.  To  take  advantage  of  the  proposed  power  query 
algorithm  it  is  therefore  necessary  to  realize  power  queries  Wf{exp{^Mq))  on 
a  quantum  computer  in  such  a  way  that  cost(Q',p)  =  o(p )  •  cost(Q/) 
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The  implementation  of  power  queries  with  cost  that  is  not  linear  in  the 
power  p  of  the  query  is  still  not  settled  and  requires  more  work.  It  would 
be  of  interest  to  identify  subclasses  Q'  C  Q  for  which  we  are  able  to  prove 
cost(Q/,p)  =  o(p)  ■  cost(Q/). 

Another  open  question  is  whether  it  is  possible  to  extend  the  methods 
we  used  for  upper  and  lower  bounds  for  the  Sturm-Liouville  problem  in  one 
dimension  to  the  Sturm-Liouville  problem  in  higher  dimensions.  Most  im¬ 
portant  for  this  problem  is  probably  the  extension  of  the  results  in  pj  on 
good  enough  approximations  of  the  eigenvector  with  the  smallest  eigenvalue 
to  higher  dimensional  Sturm-Liouville  problems. 
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A  Proofs  of  Theorem  [3]  and  Lemma  0 

Proof  of  Theorem  [3|  We  want  to  prove  that  for  a  unitary  operation  Q  with 
fixed  eigenvectors  |-0S)  and  varying  eigenvalues  e2miPs,  s  =  1  every 

quantum  algorithm  with  power  queries  Wj'(Q)  as  in  Definition  [2]  can  be 
written  as 

jV-(T)W)>  =  UTWp(Q)UT^  . . .  U,WfPQ)U^)  =  sf  V)  l*> . 

k 

In  this  equation  the  trigonometric  polynomials  S^\(p)  = 

are  of  powers  J7r,  with  Jt  recursively  in  Equation  (JTl|.  We  also  claim  that 

the  coefficients  G  C  fulfill  E/eyT  |a£j|a  =  1- 

The  proof  is  by  induction  on  the  number  of  queries  T.  We  will  write 
the  state  of  the  algorithm  after  T  steps  | in  the  basis  (| k)  \ 'f>s))k,si 
which  is  split  into  a  control  part  | k)  and  an  eigenvector  part  | ips).  We  will 
not  address  the  ancilla  qubits  in  our  proof,  but  they  can  easily  be  treated 
(after  possibly  reordering  the  qubits)  as  control  bits  that  are  never  used. 
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At  the  start  of  the  algorithm  we  can  write 


k.s 


which  contains  only  powers  elmp:p  from  j  G  Jo  =  {(0, . . . ,  0)}  and  obviously 


E  E  i«j®/  =  v 


a 


(o)  i: 

fe,s,(0,...,0)  I 

k,s  JeJo  k,s 

Let  us  now  assume  \ip^(Q)')  can  be  written  as 


=  1. 


\*mm = E  E  aT,/"10  i*>  w.> . 


with  coefficients  a\^ j  fulfilling  condition  (12).  If  we  apply  W[r'r^'  (Q)  to 


\ip(T\Q))  we  get  (kiT+1  is  the  control  bit,  i.e. ,  the  /T+1-th  bit  in  the  binary 
representation  of  k): 

Kh‘W)b(T>w)>  =  E  Eqm/",v mw 


k,s  JejT 


klT+ 1=° 


+  E  E"S/"”i‘)fl"tl I*)- 

k,s  jejT 

klT+l  =1 

We  analyze  the  second  term  where  the  control  bit  is  set  and  get  the  following 
sum  over  J7r+i  (here  e*s  is  the  s-th  unit  vector  (0, . . . ,  0, 1,  0, . . . ,  0)): 

E  <)/“E  k)Qn*'  w  =  E  d^2^2’^1’’-  i*>  w 


j&Jr 


j&Jr 

■■  a^Ke2m^+PT+ie~a)-v\k)  \4)a/ 


j&Jr 

=:  E 

jeJr+i 


(T)  p27ri(j+pT+iel)-(^ 

|fc)  |Vb) 


e27rik'v 

k,s,j 


where  is  dehned  as 

K,S,J 


a 


(T+l)  ___ 
~k,sj 


a 


a 


(T)__ 

k,sJ-pT+ iel 
(T)__ 

k,s,j 


for  fciT+1  =  1  and  j  -  pr+iel  6  Jr 

for  kiT+1  =  0  and  j  G  Jt 
otherwise 
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Thus  we  can  write 


iupxT+tw)h(T'w)>  =  E  E 

k’s  j£Jt+i 

and  our  normalization  condition  |l2j)  holds  for  a^+}\ 


EE 

k’s  j£Jr+i 


(T+l)  ,2  = 
k,sj 


E  E 

k,s  JejT 
kiT+ 1=0 

EE 


\a 


(T)j2 

k,s,j' 


E  _  E 

k,s  j-PT+ieleJT 


a 


(T)__ 

k,sj~PT+  ie1 


klT+ 1  1 


a 


cn  i2  _ 


/c,s 


k,s,j 


=  l. 


j&Jr 


The  next  step  is  to  apply  Ut+i  to  (Q) .  Define  the  coeffi¬ 

cients  uitt,k,s  =  (Ml  Ut+i  | k,s).  Then  we  can  define  a  jM 1 }  such  that 

Ut+1W%+'(Q)\^t>(Q))  =J2  E  ^‘^UT+1 1*)  |« 

jGJr+i 

-E  E  E5Ej,“-.«,.e2’i)Vii>iv-.) 

D  jGJ7r+i 

„(T+1)„27riJ-i^ 


:;E  E  rv"«  io  i*> 


1 >*  3&Jt+i 


It  remains  to  check  that 


E  E 

jeJr+1 


1  (T+l) 

a  - 

I  i,t,j 

2 =E  E 

EOOTM' 

~(T+1) 

Z_.  ak',s'JUU,k',s’ 

l’t  jeJr+i 

_ 

_k'  ,s' 

E  E 

k,s,k',s'  j£jT+1 


a 


(T+l) 

mJ 


^  '  (  iq,t,fc,s  )  U l,t,k',s ' 

L  l,t 


a 


(T+l) 
k' , s', j 


EE 

fc’s  JeJr+i 


a 


(T+l) 


=  1. 


□ 


Proof  of  Lemma [^}  Consider  quantum  queries  exp(|Mg)  for  constant  func¬ 
tions  g(x)  =  q  G  [0,  1)  in  the  Sturm-Liouville  problem.  From  equation  (13) 
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we  know  that  the  final  state  of  every  T  power  query  algorithm  can  be  written 
as 

k<T,(?)}  =  £  £ 

k  m£jVlT 

Let  B  be  a  partition  of  the  set  of  all  basis  states  |  k).  Thus  the  probability  to 
measure  a  state  from  the  set  B  G  B  is 


£ 

£ 

fees 

mGAir 

£ 

£ 

keB 

jrni&M 

£ 

£ 

(T)  -zmq 

/Ue2 


dT)  \*e -|ml9 


2e2 


(T)  „n™2  q 


kdB  mi,m2&MT 

■  £ 

1&CT 


.ni2£MT 

ST)  YJT) 


where  the  set  Ct  is 


Ct  =  {mi  —  m2  :  mi,m2  G  .Mr}  •  (25) 

and  with  coefficients  pR  /  dehned  by 

*§?  =  £  £  ('£,)*-)£,■  (26) 
/cGB  mi  ,m2 GA^t 
m2~mi=l 

(T) 

For  any  partition  £?  we  can  now  bound  the  /3R  /  as  follows 


£K?=££  £  (C,R 


B&B 


B£B  k£B  m\,m,2&MT 
m,2—m\=l 


(T) 

T  v  7 

:,m2 


< 


£  £ 

fc  m\,m2&M.T 

m,2—mi=l 


(T)  (T) 


where  is  the  sum  over  all  possible  states  h).  From  (15)  we  now  derive 
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by  the  Cauchy- Schwarz  inequality 


E  E 

k  m\,m2dMT 

m,2—mi=l 


(T)  (T) 

VlmA,™  2 


E  E 


(T)  (T) 

Vk,mVk,m+l 


<  1. 


It  remains  to  show  that  the  two  definitions  of  C+  in  equations  ([l7|)  and  ( 25 ) 
are  identical.  The  proof  is  by  induction.  T  —  0  is  trivially  true.  We  use  the 
definitions  (|ITj),  ([14])  of  JT  and  Mt  to  see  that 


Mt+  1  —  jjl  +  •  ■  ■  +  jn  '■  j  e  Jt+\  | 

=  {il  +  •  ■  •  +  in,  jl  +  ■  ■  ■  +  jn  +  PT+ 1  :  J  e  Jt} 

=  {m,  m  +  Pt+i  ■  m  G  ■ 


With  this  we  can  derive 


jCt+ i  =  {mi  -  m2  :  mi,  m2  G  Mt+ i} 

—  {mi  -  m2,mi  +  pT+ i  -  m2,mi  -  m2  -pr+i, 
mi  +  Pt+i  ~  m2  -  Pt+i  ■  mi,m2  G  Mt } 
=  {/,/+  Pt+i,  l  ~  Pt+i  '■  l  £  £t}  , 


which  completes  the  proof. 


□ 
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